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We discuss the difficulties in testing the hadron models by using the A'^* parame- 
ters extracted from the empirical amplitude analyses of the -kN and reaction 
data. As an alternative or perhaps a more advantageous approach, we present a 
Hamiltonian formulation that can relate the pion production dynamics and the 
constituent quark models of N* structure. The application of the approach in 
investigating the A and A'^*{S'ii) excitations is reviewed. It is found that while the 
A excitation can be described satisfactory, the ttN scattering in Sn channel can 
not be described by the constituent quark models based on either the one-gluon- 
exchange or one-meson-exchange mechanisms. A phenomenological quark-quark 
potential has been constructed to reproduce the 5ii amplitude. 



One of the challenging theoretical problems is to understand the hadron 
structure within Quantum Chromodynamics(QCD). There exist two different 
approaches. The most fundamental one is the Lattice QCD calculation. The 
second one is to develop various QCD-inspired models. While a lot of pro- 
gresses have been made in this direction, only very limited attention has been 
paid to the need of developing appropriate reaction theories for testing their 
predictions by using the data of ttN and reactions. In the past few years, 
we have addressed this question concerning the constituent quark models. Here 
we would like to review the progresses we have made [@-^ and discuss future 
directions. The work for N*{Sii) involves a collaboration with T. Yoshimoto 
and M. Arima. 

We were motivated by the following observation. It is common to compare 
the masses and decay widths predicted by the constituent quark modcls(or any 
existing hadron models) with the data listed by Particle Data Group(PDG). All 
calculations of decay widths have been done jJ-H] perturbatively. It has been 
found that such a perturbative calculation can at best describe the general 
qualitative trend of the data, but not the quantitative details. It is important 
to note here that the PDG's values are extracted from the empirical nN ttN 
and t:N amplitudes which contain both resonant and non-resonant com- 

ponents. In most partial waves, the non-resonant mechanisms are important ; 
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one can see this from the fact that most of the resonances identified by PDG are 
in fact not visible in ■nN and 'yN cross section data. By the unitarity condition, 
therefore the extracted resonance parameters inherently contain non-resonant 
contributions. Clearly, except in a region where the non-resonant contributions 
are negligibly small, the comparison of the PDG values (or values from other 
amplitude analyses) with the decay widths calculated perturbatively from the 
constituent quark models could be very misleading. In particular, a perturba- 
tive calculation of decay widths is obviously not valid for cases in which two 
nearby resonances in the same partial wave can couple with each other through 
their coupling with the meson-nucleon continuum. Similar precautions must 
also be taken in comparing the predicted masses with the PDG values. 

To have a more direct test of constituent quark models, it is necessary 
to develop a nonperturbative approach that takes account of the unitarity 
condition and can relate the ttN and reactions directly to the predicted 
internal quark wave functions of baryons. We have achieved this by developing 
an approach within the Hamiltonian formulation. In the following, we will 
first describe the major steps and discuss how our previous work on the 
A excitation can be interpreted within this framework. We then discuss our 
findings from an investigation of the excitation of the A'^*(5ii) resonances in 
ttN scattering. 

We start with the usual constituent quark model defined by the following 
Hamiltonian 

hB= K + V^oni + Vqq. (1) 

where K is the kinetic energy, Vconf = X]i<j otcfij is the usual linear con- 
finement potential. For the residual gg-interaction Vq^jp. Eq.(l), we consider 
both the usual one-gluon-exchange(OGE) model pO|-p^ and the recently de- 
veloped one-meson-exchange(OME) model(the weak 77-exchange is sup- 
pressed here, but it was included in our investigation ^). Explicitly, we have 

K'''' = E < • > ['^^ • ^l^'^'-'^j) + S^,V^in,)] , (2) 

K,""" = E • t*^'^ • '^:.-^/(^«.) + S^,V^i^v)] , (3) 

i<j 

where the color SU(3) factor is < • Xj >= — |, (Ti and are respectively the 
spin and isospin operators, and Sij is the usual tensor operator. The radial 
parts of the potentials in Eqs.(2)-(3) are given in Ref. We only note here 
that they are regularized by form factors F°'{q) = A^/(A^ + q^)- This is 
consistent with the notion that the constituent quarks are not point particles 
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within an effective theory. This regularization of the gq-potcntial is essential 
in obtaining convergent solutions for the bound state problem defined by the 
Hamiltonian hs (Eci-(l))- If the potentials are not regularized by form factors, 
the ground state energy is not bound from below. 

It is important to mention here that the considered two models are rather 
different mainly due to the flavor-dependent factor • Tj in OME model. This 
has important consequences in predicting the baryon spectra, as discussed 
in Ref. p^ . It is fair to say that with suitable adjustments and additional 
phenomenological parameters, both models can reproduce the general pattern 
of PDG's baryon spectra. Our objective is to find a way to distinguish them 
by considering nN and 7^'^ reactions. The situation is similar to atomic and 
nuclear physics. Only by investigating reactions, the dynamical content of the 
theoretical models can be truly tested. 

With the Hamiltonian Hb (Eq.(l)) defined above, our first step is to solve 
the three-quark bound state problem 

hB\B)^mB\B), (4) 

where \B) is the baryon wave function with the label B denoting collectively 
the spin-parity J'^ and isospin T; niB is the mass eigenvalue. We use the 
diagonalization method developed in Ref. [|l5| to solve Eq.(4) by expanding 
the baryon wave function as 

|i3)_|j-T)=^af^|J-r;z), (5) 

i 

where the basis states are appropriately antisymmetrized harmonic oscillator 
wavefunctions. The coefficients ^ in Eq.(5) and the mass eigenvalues mj^T 
of Eq. (4) are obtained from diagonalizing the matrix 

m,, = {J^T-i\hB\J^T-j). (6) 

In practice diagonalization is performed within a limited number of basis states. 
Then the solution of Eq.(4) is a function of the oscillator range parameter b. 
We treat it as a variational parameter and find h by imposing the condition: 

The basis state is chosen so that the mass eigenvalue mjirj' does not change 
by further extension of the basis states. In practice we include the basis states 
up to ll?ia;. 
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The next step is to define how an external meson(M) can interact with the 
bound three-quark systems. Namely, we need to calculate the matrix elements 
for the transitions from one-baryon states(i3) to meson-baryon states(Af'i3') 



{B'-M'\HM{k)\B) 



(8) 



where B is a bound state wave function generated from the above structure 
calculation, and HM{k) is an appropriate operator describing how a meson M 
with a momentum k is emitted or absorbed by constituent quarks. Following 
the previous works we assume that HM{k) for M = tt is a one-body 

operator which can be derived directly from the nonrelativistic limit of the 
Feynman amplitude ~ ^^^upij^j^k^j^up for the q irq transition. To be 
consistent with the nonrelativistic treatment of constituent quarks, we keep 
only the terms up to the order of p/mq. In coordinate space, the resulting 
q + TT q transition operator is 



F{k), (9) 



where a denotes the z-component of pion isospin and (p^) is the derivative 
operator acting on the initial (final) baryon wave function; k and uJt: are the 
momentum and energy of pion, respectively. 

We fix the irqq coupling constant by assuming that the ttNN vertex func- 
tion can be calculated from Eqs.(8)-(9) using the (Os)^ nucleon wavefunction. 
We then find that — | , where Ad^ is the observed mass of the nu- 
cleon and we use the empirical value g^j^^/Air — 14. The same /jr^g coupling 
constant is also used to evaluate the qq potential Eq.(3) of OME model. Ac- 
cordingly, we also introduce a form factor F{k) in Eq.(9) to account for the 
effect due to the finite size of constituent quarks. 

For photon coupling, we calculate the B B'j by assuming the commonly 
used J — qq vertex interaction. The resulting form is 



{k)^{B';^\H^{k)\B), 



(10) 



with 



H-,{k) 



(11) 



By using the (Os)^ wave functions for N and A and setting ^ 300 MeV 
to evaluate Eqs.(10)-(11), one finds that the nucleon magnetic moments can 
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be reproduced well and the resulting helicity amplitudes for the A 

transition are A3/2 160 x lO'^ {GeV)-^/^Ai/2 = A3/2/V3). These values 

of helicity amplitudes, which are about 40 % lower than the values listed by 
PDG but are close to various quark model predictions are consistent 

with our interpretation in Ref. that the constituent quark model predictions 
can only be compared with the bare values extracted from the data by removing 
the on- and off-shell effects due to non-resonant interactions in a dynamical 
approach based on the Hamiltonian formulation. 

With the vertex functions defined above, we can then define a hadronic 
Hamiltonian for investigating ttN and 'yN reactions. It has the following form 

H = Ho + Hj, (12) 

where the free Hamiltonian can be written in a second-quantization form 

Ho = Y, f dpeB{p)bl{p)bB{p)+J2 I dkujM{k)a\.j{k)aM{k). (13) 

Here, &^(&s) and aj^^(aM) are the creation (annihilation) operators for the 
baryons and mesons respectively, and ebIp) = {m% +P^Y^'^ ,^M{k) = {m\j + 
fc^)^/^ . It is important to note here that the baryon mass tub in Eq.(13) 
is generated dynamically from solving the three-quark bound state problem 
defined by Eqs.(4)-(7). Their values could be significantly different from reso- 
nance positions listed by PDG. We use the experimental value for the meson 
mass toa/. 

The interaction term in Eq. (p^ ) is written in terms of the vertex functions 
defined by Eqs.(8) and (10) 

Hi=Y. [ dpdp'dk[{B'\HM\B-M)b^B'iP')bBiP)aM{k)+h.c^. (14) 

BB'M 

The above interaction Hamiltonian is similar to that of the dynamical model 
developed in Ref. [|l|, except that the B' EM transition amplitudes are 
now determined by the predicted quark wavefunctions. As discussed in Ref. 
10, it is a non-trivial many-body problem to calculate nN and reactions 
with the use of Hj. To obtain a manageable reaction theory, we follow Refs. 
|^,[l6|] and apply the unitary transformation up to the second order in Hj to 
derive an effective Hamiltonian. The essence of the unitary transformation 
method applied in Ref. ^ is to absorb the unphysical transition B M' B' 
with TUB < mB> + TTiM' into non- resonant potentials. The resulting effective 
Hamiltonian then takes the following form 

H,s = Ho + r + T^ + V, (15) 
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where Hq is defined in Eq. (^3|). The vertex contains only the physical 
decay process B — > M' B' with niB > mM' + "ms' 

E j dkdpdp' {B'-M'\HM,,\B)h^B'{P>\vik')bB{p) 

MBB' 



X 6{mB - {ruB' + mM')), (16) 

where e{x) = 1(0) for a; > 0(a; < 0). The non-resonant MB M'B' two-body 
interactions are defined by 

I dkdk'dpdp\B';M'\v\B-,M)a\j,{k')aM{k)b\,{p')bB{p).{ll) 

MM'BB' •' 

The construction of the matrix elements of the nonresonant potential v is 
explained in Ref . ^ . Typically, it consists of the energy-forbidden s-channel 
terms such as ttN ^ TV — > ttN and vrA ^ A ^ ttA, and particle-exchange 
terms such as the usual crossed nucleon and A terms and the p- and w-exchange 
terms. 

By using the standard projection operator method 0, it is straightforward 
to derive from the effective Hamiltonian Eq. ( |l5| ) a calculational framework for 
ttN and reactions leading to various final meson-nucleon final states. The 
resulting transition operator can be written as 

T^,p = tc.,p+Y.'^lN: [D-\E)]^^^tM;,p. (18) 

Here a,/? denote the meson-baryon states such as ^N, nN,rjN and ttA. N* 
are mass eigenstates of Eq. (^). The first term in Eq. (|l^) is the non-resonant 
amplitude determined only by the non-resonant interaction v 

ta,f3 ^ Va,p +Yva,-tG"{E)t^i3, (19) 
7 

with 

[G'^{E)] -'^E-sbM- ^M-,{k) + ze. (20) 



The second term in Eq. (|l^) is the resonant term determined by the dressed 
A^* propagator and the dressed vertex functions: 

[D{E)]^^^ = {E- mN^)6,, ~ E,,,(i?), (21) 

f AT'.a = ^ r^*.7 ['^7" + ^7(^)^7, a] J (22) 

7 

f U; = E ['^70 + t^,',G°{E)] tI^,. (23) 
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Figure 1: The recent Mainz data[19] of differential cross section and photon asymmetry S 
of the p(7, ■iT^)p reaction are compared with the predictions from the SL model. Ge is the 
strength of the E-^^-^- amplitude of the A 7A'' transition. 



In Eq. (21), the A^* self-energy is defined by 



^^AE) - Vr^;,7G°(£;)ft (24) 



We now note that the above coupled equations relate the full scattering 
amplitude Ta^p Eq.(18) nonperturbatively to the quark wavefunctions calcu- 
lated from the considered constituent quark model. In particular, the decay 
widths, as hsted by PDG, correspond to the dressed vertex TN*,-y defined by 
Eq.(22), and hence they are expected to be different dynamically from the 
bare vertex functions Tn-^^ which are calculated from using Eqs.(8)-(11). It 
is also clear that the accuracy of the non-resonant interaction Va,[3^ which de- 
termine the nonresonant amplitude ta.p via Eq.(19), plays an important role 
in identifying the resonant amplitude of Eq.(18) from the data. It is therefore 
also essential to have it calculated from the same quark model. This can be 
achieved in our approach through the use of unitary transformation method, 
as discussed above(the developed method can be straightforwardly extended to 
account for the interactions due to heavy mesons). Without such a consistent 
treatment of both the resonant and non-resonant interactions, it is difficult to 
draw conclusions. 

We now turn to discussing our findings. First we consider the A state 
in P33 channel. For this investigation, we only have one N*{= A) state and 
two channels a,(3 = ttN,jN. Eqs.(18)-(24) are then reduced to the set of 
equations developed in Ref. (called SL model from nowon). The accuracy 
of the SL model are illustrated in Figs. 1-2. We see that the predictions are 
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Figure 2: The MIT-Bates data[20] of asymmetry Ttl of the p(e, eVO) reaction at = o.l26 
(GeV/c)^ are compared with the predictions by the SL model. 



in good agreement with the recent data from Mainz and MIT-Bates pOt . 
It is therefore sufficient to adjust the parameters characterizing the the qq 
potentials Eqs.(3) or (4) and the constituent quark form factor F{k) for the 
vertex interaction Eq.(9) to fit the bare A parameters of the SL model. This is 
a rather restricted procedure since there are only about 4 parameters for each 
of the considered OGE and OME models. 

In Fig. 3, we show that the A — > ttN form factor calculated from using 
Eqs.(8)-(9) agree to a large extent with the phenomenological bare form fac- 
tor (dashed curve) of SL model. The predicted mass (via Eq.(4)) is also required 
to be the bare mass toa = 1300 MeV of SL model. 

Once the parameters of the considered constituent quark model are de- 
termined, we then can generate the mass m,B and quark wavefunctions for all 
possible partial waves by solving Eqs.(4)-(7). The relevant vertex functions 
can then be calculated by using Eqs.(8)-(11) and the TriV and reaction am- 
plitudes can be predicted by solving Eqs.(18)-(24). The comparison with the 
data will tell us whether the considered constituent quark model is correct. We 
have explored this by investigating the TriV scattering in partial wave. The 
calculation was done by considering two N* and three channels ttN, rjN and 
IT A. The predicted ttN phase shifts are shown in Fig. 4a. We see that neither 
OGE(solid curve) or OME(dashed curve) models can describe the data(open 
circles) in the entire energy region. 

As an attempt to improve the fit to the 6*11 amplitude, we have also ex- 
plored the mixture of OGE and OME models. It turns out that such a hybrid 
model also fails, mainly due to the very disruptive tensor component of the 
OME model in determining the phases of wave functions. 
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Figure 3: The A — > vrA'^ vertex function(solid curve) determined within the OGE model is 
compared with the phenomenological form(dashed curve) of the SL model[l]. 

We therefore turn to investigating a purely phenomenological model. By 
analyzing an analytical model |^ , we have found that the data of A excitation 
and Sii ttN scattering seem to favor a tensor term due to one-gluon-exchange 
and a spin-spin interaction due to one-meson-exchange. This has guided us to 
explore many phenomenological models. For example, we have found that the 
ttN Sii amplitudes can be much better described by the following phenomeno- 
logical model 



The results are shown in Fig. 4b. We see that the general feature of the data 
can now be reproduced. The remaining discrepancies perhaps can be removed 
if we refine various form factors in the qq potential and N* — > nN, rjN, tt A 
vertex functions(Eq.(9)). 

In conclusion, we have developed a Hamiltonian formulation for relating 
the dynamics of TriV and 7iV reactions to the constituent quark models. The 
approch has been applied to investigate the A and iV*(S'ii) resonances. It is 
found that while the A excitation can be described satisfactory, the TriV scat- 
tering in Sii channel can not be described by the constituent quark models 
based on either the one-gluon-exchange(OGE) or one-meson-exchange(OME) 
mechanisms. The data seem to favor the spin-spin interaction due to one- 
meson-exchange and the tensor interaction due to one-gluon-exchange. A phe- 
nomenological quark-quark potential has been constructed to reproduce the 
Sii amplitude. 



(25) 
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Figure 4: The predicted -nN phase shifts 
from OGE(sohd curve) and OME(dashed cu 
Eq.(25). 
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are compared with the data(open circles): (a) 
rvc) models, (b) from phcnomcnological model 



To end, we would like to point out that the multi-channels multi-resonances 
parameterization of TriV reaction amplitudes, first developed by Cutkosky and 
his collaborators and has been revived recently [|l8|, can be derived pH 
from Eqs.(18)-(24). The bare parameters associated with this phenomenolog- 
ical model can have dynamical interpretations within our Hamiltonian formu- 
lation. Consequently, it may be more advantageous to apply this parameter- 
ization to analyze the forthcoming data. The extracted bare iV* parameters 
can be used as the data for exploring the constituent quark interactions using 
the procedures detailed in Ref. |^ and outlined in this contribution. 
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